Exploring the origin of the fine structures in the CMB temperature angular power 

spectrum 
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The angular power spectrum of the cosmic microwave background (CMB) temperature 
anisotropics is a good probe to look into the primordial density fluctuations at large scales in the 
universe. Here we re-examine the angular power spectrum of the Wilkinson Microwave Anisotropy 
Probe data, paying particular attention to the fine structures (oscillations) at t = 100 ~ 150 re- 
ported by several authors. Using Monte-Carlo simulations, we confirm that the gap from the simple 
power law spectrum is a rare event, about 2.5-3cr, if these fine structures are generated by exper- 
imental noise and the cosmic variance. Next, in order to investigate the origin of the structures, 
we examine frequency and direction dependencies of the fine structures by dividing the observed 
QUV frequency maps into four sky regions. We find that the structures around I ~ 120 do not 
have significant dependences either on frequencies or directions. For the structure around i ~ 140, 
however, we find that the characteristic signature found in the all sky power spectrum is attributed 
to the anomaly only in the South East region. 



I. INTRODUCTION 

The inflationary cosmology is a successful paradigm 
in explaining the generation of primordial density fluc- 
tuations and solving essential problems of the classical 
Big Bang cosmology The primordial density fluc- 

tuations are transformed into the anisotropy of the Cos- 
mic Microwave Background (CMB) and the Large Scale 
Structure of the universe (LSS). Thanks to the high 
angular resolution and longer-term observations of the 
anisotropy of the temperature fluctuations, such as by 
the Wilkinson Microwave Anisotropy Probe (WMAP) 
0, the South Pole Telescope [|[, the Atacama Cosmol- 
ogy Telescope M , the Arcminute Cosmology Bolometer 
Array Receiver [10| . the Cosmic Background Imager [ll[ 
and so on, it has been found that the angular power spec- 
trum of temperature fluctuations conforms to the predic- 
tion from the ACDM model with slow roll inflation. 

While the observed angular power spectrum is globally 
consistent with a smooth power-law primordial spectrum 
of density fluctuations [l3 ? - fl6j , some gaps between the 
prediction from the simplest power-law model and the 
observed data have been reported. They are discussed 
recently by grace of the detailed observations, including 
a small bump and dip at t = 20 - 40 [l7Hl9| or oscillation 
around I = 100 - 150 [Mil- 

Possible causes of these anomalous structures are dis- 
cussed in the literature. The small bump and dip at 
£ = 20 - 40, for example, may be explained by the mass 
variation of the inflaton during the inflation phase (l7l - 
fl9L l24j . When inflaton obeys single slow- roll inflation dy- 
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namics, the power spectrum of primordial density fluctu- 
ations shows a power-law feature, and the angular power 
spectrum of the CMB is expected to be a smooth curve. 
If the inflaton mass has changed during inflation, how- 
ever, some oscillating structures emerge in the power-law 
primordial power spectrum because the inflaton field is 
forced to accelerate and/or decelerate rapidly during the 
slow roll inflation phase. In such cases the bump and dip 
structure rises up in the angular power spectrum of the 
CMB. 

On the other hand, it seems difficult to explain the 
oscillating structures around I = 100 - 150 on the firm 
theoretical background, though some works have tried to 
explain the origin [20H23l |. In the previous paper, some 
of us have tried to explain the oscillating structures with 
the inflaton mass variation [23}. The condition consid- 
ered in that paper is that inflaton mass changes with 
oscillations. We found that oscillating structures can be 
generated at the arbitrary scale by adjusting oscillation 
number and time scale of the mass variation. However, 
the width of this structure tends to become so wide, and 
to match the observed data we need some fine-tunings. 
If we force to explain the oscillation matching with the 
observed data, the parameters become unrealistic values. 
Therefore, we have concluded that it is difficult to explain 
the oscillating structure on the angular power spectrum 
found at multipole range I = 100-150 with the inflaton 
mass variation. 

Nakashima et al. [22j have proposed a sudden change 
of the sound velocity of the inflaton field during infla- 
tion. Based on their model, the authors found oscillat- 
ing structures in the primordial power spectrum. How- 
ever, the oscillating structures tend to extend in a wide 
range of wavenumber up to I < 300 which includes the 
first acoustic peak. Therefore, this model may not match 
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with observed data which shows the oscillation only at 
the confined region of 100 < I < 150. 

In light of the difficulty in explaining the structure on 
the theoretical background, in this paper we closely ex- 
plore the origin of the structure with the data taken by 
WMAP, particularly paying attention to the frequency 
and direction dependences. If the observed structure is 
really from the cosmological origin, we expect the struc- 
tures should be independent of them. 

This paper is organized as follows. In section II, we 
review the method to estimate the angular power spec- 
trum from the two point correlation function. Following 
the method, we examine the probability of the fine struc- 
tures being generated from noise and/or cosmic variance 
in section III, and estimate the angular power spectrum 
with each frequency band in section IV. In section V, we 
give the angular power spectrum with the partial sky, 
with an explanation of some difficulties in estima ting the 
angular power spectrum with the partial sky map [25l . l26j . 
Section VI is devoted to the summary of this work. 



II. THE ANGULAR POWER SPECTRUM 

The angular power spectrum Ce is a good indicator for 
quantifying the temperature fluctuations of the CMB, 
and defined as 



(a* lm ai> m i ) — Su'Smm'Ce 



(1) 



with 



#AT(M)r/ m (0,0), (2) 



where Y^ m (6, (ft) is the spherical harmonics function eval- 
uated at the position (8, eft) of spherical coordinate, 
AT (6, (ft) is the value of temperature fluctuations and 
(...) means ensemble average. In practice, however, we 
can observe only one unique sky and need to estimate 
the angular power spectrum from one realization. We 
can obtain the estimator Ce of the Ce by 



C, = 
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(3) 



In this work, we use a fast CMB analysis named Spice 
introduced by Szapudi et al. for temperature (27j and 
Chon et al. for polarization [28|. This software calcu- 
lates the angular power spectrum from the CMB tem- 
perature/polarization data and can convert the angular 
power spectrum to two point correlation function or vice 
versa, if we need. We review the relation between the an- 
gular power spectrum and the angular correlation func- 
tion for later discussion. 

The CMB temperature angular correlation function 
£(0) is defined by 



where n and n' denote sky directions, and cos 6 = n ■ 
v! . The relation between the angular correlation function 
and the angular power spectrum can be written as, 



(5) 



Ci = 2tt sm6d9Z(d)P e (cosd) , (6) 



where Pe(cos8) is the Legendre polynomial. 



III. VALIDITY OF THE FINE STRUCTURES 
A. Monte-Carlo method 

As a possible origin of the fine structures, instrumen- 
tal noise and cosmic variance effects should be called 
into question first. In this section, we examine the pos- 
sibility of the fine structures being generated by these 
noise effects. For this purpose, we generate a number 
of sky maps of CMB temperature fluctuations using the 
HEALPix function, isynfast [29| . The routine generates 
a temperature fluctuation map from an input angular 
power spectrum Ce- For the Ce, we adopt the best-fitting 
ACDM model of the WMAP seven year parameter table 
[30| . The resolution of this map is taken as N s id c = 512. 
This value is same as that of the released sky map given 
by the WMAP team. 

Next, we add the Gaussian noise expected for the 
WMAP observation on these simulated sky maps. The 
variance of the instrumental noise can be written by 



(7) 



bs 



ae) = (AT(n)AT(n')] 



(4) 



where superscript i represents the pixel number on the 
sky, cr the rms noise per an observation, and N* hs the 
number of observation of zth pixel. The values of do and 
iV* bs are also given by the WMAP team Q . In this way, 
we can prepare the sky maps of CMB temperature fluctu- 
ations of different universes with the same cosmological 
parameters and noise properties. 

In the present analysis, we prepare 3,000 sky maps and 
calculate the angular power spectrum for each sky map. 
These Cg's are denoted by C|'s using a superscript of the 
map number j . We estimate the average and variance of 
these C^'s at each multipole moment t, and compare the 
angular power spectra with the WMAP angular power 
spectrum, cW MAP . j n orc { er ^ emphasize the differ- 
ences from the ACDM model, we show the distributions 
of the simulated Cg's as residuals from the ACDM model, 
C-ACDM in Fig _[jj The variances f q> s at 1(7j 2c r, and 3a 

are shown as the boxes. We can see that the average of 
the simulated Cfts (the black line in the figure) is on the 
zero horizontal axis. The red line represents the residuals 
between C £ ACDM and C £ WMAP . 
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FIG. 1. Differences between the WMAP angular power spec- 
trum C| WMAP and the power spectrum of the ACDM model 
C* CDM (red line). The boxes represent the variance of 1, 2, 3a 
at each I estimated from the 3,0000 simulations. The black 
line represents the average value of residuals between Cj 's 



and C, 
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FIG. 2. Contours of equal probabilities as a function of N ltT -2 a 
and N2<y. The red, green and blue lines correspond to la, 
2a and 3a, respectively. The (yellow) point represents the 
WMAP data. 



B. significance of the structure 

In the multipole region of £ = 100 - 150, the distribu- 
tions of the CYs are well approximated by the Gaussian 
one. In that case the probability that the Ce is within 
la and 2a should correspond to about 68.3% and 95.5%, 
respectively In our analysis, there are 51 data points 
because we focus only between £ = 100 - 150, and the 
expected number of data over la or 2a is estimated as 
16.2 or 2.3 points, respectively. Thus, it is natural that 
some data points deviate from the theoretical curve to 
this extent. However, from Fig. [TJ we find that 19 and 
7 points of data are over la and 2a levels, respectively. 
That indicates the Cj^ MAP somewhat deviates from the 
standard Gaussian distribution more than expected. 

We roughly estimate the probability P as a function of 
N 2 a and Nx a - 2a which are the numbers of realization that 
exceeds 2a and lies between la and 2a, respectively, as- 
suming the realization follows the Gaussian distribution. 
We show the result as contours in Figure. [2J In the fig- 
ure the red, green and blue lines represent la, 2a and 3er, 
respectively. From the figure, we can see that the peak 
is at the expected value (Ni a - 2cr , N 2a ) = (13.9, 2.3), and 
the probability decreases with distance away from the 
expected value. The yellow point represents the value 
from the WMAP data; (N la - 2cr , N 2a ) = (12, 7). It lies 
at the location between 2.5-3cr, which is consistent with 
the result of [2(| ■ Therefore, we can understand that the 
fine structures at the range of £ = 100-150 of the angu- 
lar power spectrum which are observed by WMAP team 
seem a rare event, if these structures are generated only 
by the noise and the cosmic variance effects. 

If the structures were not related to the noise nor cos- 
mic variance, the other possibilities would be foreground 



effects and/or real cosmological signal. From the next 
section, we examine frequency and direction dependences 
of the fine structures in order to investigate whether the 
structures come from cosmological or astronomical phe- 
nomena. 



IV. FREQUENCY DEPENDENCE 

If any astronomical phenomena such as syncrotron 
emission, dust emission and radio galaxies infiltrate the 
CMB temperature fluctuation data and create the fine 
structures, we expect that they have a characteristic fre- 
quency dependence. Therefore, we estimate the angular 
power spectra for three different frequency bands, namely 
Q, V and W bands. We find that the shapes of the fine 
structures at each frequency band are very similar to each 
other and to the all sky one (Fig. [T]) . Thus it is improb- 
able that the fine structures originate from astrophysical 
phenomena nor objects, because they will have different 
frequency dependences from the CMB blackbody spec- 
trum. 



V. ANALYSIS WITH THE PARTIAL SKIES 

In this section, we look into the differences in the an- 
gular power spectra for different sky directions. For this 
purpose we prepare four masks, namely the North West 
mask (0 < <j> < n, < 6 < tt/2), the North East mask 
(tt < </> < 2tt, < 9 < tt/2), the South West mask 
(0 < < 7T, n/2 < 9 < 7r) and the South East mask 
(tt < if) < 2ir, tt/2 < 9 < tt). Wc show these masks in 
Fig. El 
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the product to go to zero smoothly [27|, |28[. Angular 
power spectrum is then given by 



FIG. 3. Masks of the partial skies used in our analysis in the 
galactic coordinate. The left top panel refers to North West 
mask, the right top panel to North East mask, the left bottom 
panel to South West mask and, finally, the right bottom to 
South East mask. 



Using the masks, first let us calculate the gap between 
the angular power spectra of the simulated maps and that 
of the observed map by WMAP for the partial skies. Fur- 
ther, we compare the partial sky maps with each other, in 
order to see the dependence of the angular power spec- 
trum on the sky directions. However, some difficulties 
arise when analyzing the CMB map with a partial sky 
mask. Because a mask narrows down the effective area 
of the spherical surface, effective information is reduced. 
Furthermore, if a mask has a sharp cutoff at the bound- 
ary, it will lead a truncated correlation function. Com- 
puting the angular power spectrum from the truncated 
correlation function causes spurious oscillations inherent 
to the Fourier transformation, which is known as the 
"Gibbs phenomenon". However, these oscillations have 
nothing to do with the anisotropy of the CMB nor the 
fine structures, and they can be removed by apodization 
and binning techniques as we describe below. 

In the next two subsections, we review these techniques 
that we used in order to lighten and avoid these nuisance 
effects, and in the final subsection, we show the direction 
dependence of the CMB angular power spectrum. 



A. Apodization 

In fact, we find that spurious oscillatory structure in 
the angular power spectrum that arises when the masks 
are used is caused by the large amplitude oscillations in 
the two point correlation function at large angular scales. 
The oscillations are caused by the statistical errors. In- 
deed, if we use a map with a maximal angular size # max , 
nothing can be known about the correlation function for 
9 > #max- However, simply truncating the correlation 
function at 9 = 9 max does not help the situation as we 
mentioned above. Instead, the technique uses an appro- 
priate function, called the "apodization function" F(9), 
and multiplies it with the correlation function £(0) for 



C 



apd 



2tt 



0=0 



sin 9d9 £,(9)F{9)P e (cos 6) , (8) 



where # max is the maximum angle set by the mask. We 
can adopt any function as an apodization function which 
takes F(9) = 1 at 9 = and decreases as 9 increases. The 
Gaussian type Fq(9) and Cosine type Fc(9) are often 
adopted, and they are defined as 



F G (9) =cxp 



2(V81n2(T i pd) 2 



F C (9) 



COS 7T 



Capd 



(9) 



(10) 



where <r ap d = 7r0 ap d/18O° and ap( j represents the angle 
in degree where the apodization function becomes close 
to zero. Typical value of # ap d should be close to the cut 
off angle 6> max . 

We calculate the angular power spectrum with 
apodization and show the results in Fig. 2) In the upper 
panel of the figure, in the case of 9 max = 100° when 
the North West mask is used, the oscillatory feature 
with a large amplitude can be seen on the angular power 
spectrum, because we truncate the integral in Eq.(7) at 
9 = 100°. In the case of 9 = 180° when KQ75 mask is 
used, on the other hand, the oscillation is suppressed be- 
cause £ decreases as 8 — > 180°. The lower panel shows the 
results with apodization using the Gaussian and Cosine 
type functions, setting 6> ap d to 100°. As is shown in the 
lower panel, the resultant power spectra become smooth 
compared to the case without apodization (see the red 
line in the upper panel). Also it is noted that the shape 
of the apodized spectra is similar to the power spectrum 
of the full sky (6* max = 180°; the black line), though small 
scale oscillations have smoothed out. This is a bad news 
because we are interested in the fine structures. 

As we take # ap d smaller, more significant the smooth- 
ing effect for the fine structures in the angular power 
spectrum becomes. Therefore we should take the angle 
6* ap d as large as possible to keep the fine structures. The 
maximum apodization angles # ap d for our masks with 
a sufficient suppression of the spurious oscillations are 
found to be 140° and 100° for the Cosine and Gaussian 
types, respectively. In Fig. [5] we depict the apodization 
functions with those # ap d- From the figure, we can notice 
that the Gaussian type reduces the information in the 
two-point correlation function at large angles more than 
the Cosine one. In this work, we should keep the infor- 
mation as much as possible, because we focus on the fine 
structures. Therefore, we use the Cosine type apodiza- 
tion function in the following analysis. In this case, the 
structure at £ = 100-150 can survive the apodization be- 
cause the amplitude of the apodization function at this 
scale is over 0.99. 
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FIG. 5. The apodization functions F(cj>) of Cosine type (solid) 
and Gaussian type (dashed). The # ap d parameters are fixed as 
they suppress enough the spurious oscillations (#apd — 140°, 
and #° p a d SS = 100°). These lines indicate that the information 
lost is heavier for the case of Gaussian type apodization than 
that of Cosine type. 
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FIG. 4. Angular power spectra estimated from two-point cor- 
relation functions. In the upper panel, the red line represents 
the power with the North West mask and with # max = 100° 
in Eq.©, (which corresponds to a truncated correlation func- 
tion) and the black line the case with KQ75 mask. The theo- 
retical curve of the ACDM model is also shown (the magenta 
line). When we naively analyze the partial sky, the large spu- 
rious oscillations arise on the angular power spectrum. In 
the lower panel, we show the cases with apodization, with 
#a P d = 100° . The black and magenta lines are same as in the 
upper panel. The red and blue lines represent the apodized Ce 
with Gaussian and Cosine types, respectively. For both cases, 
the North West mask is used. The suppression of the spurious 
oscillations is clearly seen due to the apodization function. 



As the loss of information about high frequency struc- 
ture in the angular power spectrum is significant, com- 
ponents of Ci become correlated with each other. This 
effect and coping technique are discussed in the next sub- 
section. 



B. Binned angular power spectrum 

In the previous subsection, we see the apodization 
technique suppresses the high frequency oscillations. 
However, this technique is unavoidably accompanied 
with the loss of information. The loss appears as cor- 
relations between two different multipole moments. In 



order to see the degree of the correlations, let us calcu- 
late the covariance matrix Cu> , 

1 1 N 

C "' = c^n £ (°2 ~ SO {c\, - cv) , (11) 

where 

Q = l£q. (12) 

Here, C\ is the estimated angular power spectrum from 
the ith Monte Carlo simulation sky, N is the total number 
of the samples. From this matrix, we can estimate the 
strength of correlations between each multipole moment. 
If there are no correlations, which means the Cg can be 
estimated independently, the covariance matrix should 
be diagonal. 

On the contrary, the covariance matrix has off-diagonal 
elements if the Ci depends on another multipole moment 
component Cy. The correlations between the CYs are 
caused by the mask. If a mask is applied on the CMB sky, 
the estimated angular power spectrum is given by a con- 
volution between the power spectra of the mask and the 
true temperature anisotropics that we want to estimate. 
The convolution generates correlations between the mul- 
tipoles of A£ ~ 2-3, which makes it complicated to es- 
timate the statistical significance of the angular power 
spectrum. A simple way to obtain independent observ- 
able is to bin the angular power spectrum with the com- 
parable bin width. 

The covariance matrices of the angular power spec- 
trum with KQ75 mask arc depicted in Fig. [BJ The up- 
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per panel in the figure shows the matrix for the case 
without apodization. Each component of this matrix 
is practically vanishing except for the diagonal ones. 
This result indicates that the correlations caused by the 
KQ75 mask are not significant at the multipole region 
of 100 < £ < 150. This is because the condition that 
A£ > ~ 1 can still be satisfied with the mask, where 
$max is the maximum separation angle in the pixel do- 
main. On the other hand, in the middle panel, we use the 
KQ75 mask with apodization. In this case, the compo- 
nents around diagonal elements do not vanish. This man- 
ifests the information loss due to apodization, though the 
variances (diagonal elements) become smaller values. 

The correlations appear only between the neighboring 
multipoles. To reduce the correlations we should gather 
the neighboring components in one component, namely, 
bin the data. The contributions from the neighboring 
modes £' to the pivot scale £ are about 50% (£! = £ ± 
1) and -10% (£' = £ ± 2,1 ± 3), which is shown in the 
middle panel of Fig. |6l We find that a binning with 
A£ = 3 is sufficient because the positive and negative 
correlations nicely compensate with each other to reduce 
the correlations between the binned data. In this case 
the number of data points reduces to 17 from 51, and 
this reduction could lead some information loss. The 
covariance matrix of this case is depicted in the lower 
panel of Fig.[BJ Looking at this figure, we can see that the 
correlation becomes weak enough thanks to the binning. 
Then, each components can be considered approximately 
independent. 

Although binning is useful to simplify the statistics, 
some information should be lost in the process, especially 
if there exist fine structures in the data. In the next 
subsection we show how the binning procedure affects the 
fine structures in multipole £ = 100-150 we are interested 
in. 

We calculate the binned Ce 's and verify the fine struc- 
tures using the same steps we followed in Sec. IIII1 There 
are three ways to select the pivot scale for binning. 
We calculate the binned spectra with the three patterns 
c f = Ei=_i CW3 for case I (£ = 99 + 3ra), case I 
I (£ = 100 + in), and case HI (£ = 101 + 3n) where 
n = 1, 2, ■ ■ • , 17. The results are shown in Fig. [Jj The 
right panels represent the probability that the fine struc- 
tures come from the noise and the cosmic variance effects, 
as we have shown in Fig. [5J The result slightly depends 
on the binning pattern, but in any cases the significance 
is above 2.5 a. In the left panels of Fig. [7J we show 
the residuals between binned Q WMAP and C £ ACDM as red 
lines. The black line is the average of Ce with the 3,000 
simulations. The dispersion of the simulation data is also 
shown as the grey boxes. We can still see the oscillatory 
feature that is found in the case without binning for all 
sky, and some points are over la box or 2a box. The 
numbers of data points above 1 or 2 a are larger than 
expected from Gaussian distribution for all cases. These 
results indicate that the fine structures can survive even 
if we apodizc and bin the angular power spectrum. 
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FIG. 6. Covariance matrices of angular power spectra Cu 1 ■ 
The upper and middle panels show the covariance matrices 
without apodization and with Cosine type apodization, re- 
spectively. Non-zero off-diagonal components indicate the 
presence of correlation between I and £' . In the bottom panel 
we show the covariance matrix of the binned angular power 
spectrum. We can see that the correlation between compo- 
nents becomes weak. 
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We consider these three power spectra as the fiducial 
binned power spectra of WMAP for the full sky. To figure 
out the origin of the fine structures we compare them 
with the spectra from the partial skies obtained with the 
same apodization function and binning, as we shall show 
in the next section. 



C. The angular power spectrum with partial sky 

In the previous sections we have seen that the apodiza- 
tion and binning can successfully remove the spurious os- 
cillations and the correlations between multipoles due to 
the partial sky mask. Here we apply the same technique 
to the four different partial skies and see whether there 
exists the directional dependence of the fine structures 
in the angular power spectrum. In the following discus- 
sion, we separate the structures into two characteristic 
structures. First one is the oscillatory structure around 
£ = 100 - 120. The two peaks have a significant devia- 
tion from the smooth spectrum at around 3cr. The other 
is the bump around £ = 140. 

We show the results for the binning case I in the up- 
per panel of Fig. [3] The primary effect of the partial sky 
masks can be seen in the error bars. The estimated er- 
ror bars depend on the number of available modes TV as 
l/VN. In our analysis the effective sky coverage for the 
partial sky becomes about one-quarter compared to the 
all sky map. Therefore the error bar should be twice of 
the case of the all sky (Fig. [7|) and we can confirm this 
fact from the figure. 

We find a distinct anisotropic structure which is above 
3a error bar at £ = 138 in the South East area (bottom 
right panel of Fig. |8]), which is highlighted by a dou- 
ble circle in figure. Also, the histogram of the difference 
between ACDM and the simulated data set at £ = 138 
is shown in the lower panel of Fig. [8] for each direction. 
The vertical red lines show the differences of the power 
between the ACDM and the observed one by WMAP. 
The anisotropic structure in the South East can be seen 
in the other binning cases, as in Figs.O andfTUl 

The peculiarity of South East area as a whole at the 
multipolc range of 100 < £ < 150 can be quantified with 
the value of x 2 , defined by 

X 2 = E (Of - Cf) Cf, 1 (Cf - Cf) /Of Of , (13) 

where Cf e , is the inverse of the covariance matrix of the 
binned angular power spectrum. We show the probability 
distribution function of \ 2 , written by 

nX 2 ;n = 17) = ^ m e^/\ (14) 

and the values of x 2 for the four partial skies in Fig. [TT] 
The value of x 2 away from the peak position indicates 



the overall deviation of the observed angular power spec- 
trum from the average value of the mock power spectrum. 
From this figure, the South East area has especially pe- 
culiar x 2 value. The probability that the x 2 takes larger 
value than the South East area by chance can be esti- 
mated as 1.50 (case I), 1.56 (case II), 8.70 (case HI) %. 

These results may suggest a possibility of existence 
of characteristic features only at South East area. The 
structure which comes from cosmological origin can be 
assumed to be isotropic, and therefore this anisotropic 
structure around £ ~ 138 might be attributed to some 
astronomical origin, which have the scale about 0.6°. It 
may be interesting to note that there has been a report 
about a power spectrum anomaly around the third acous- 
tic peak at the same sky direction [2(|, although the au- 
thors argued that the origin could be the WMAP instru- 
mental noise because the third acoustic peak is located 
at the limit of the WMAP angular resolution. 

The oscillatory feature (peak and dip) that is found in 
the all sky analysis around £ = 100 - 120, on the other 
hand, can be found clearly in the North West area, and 
also found at all the other directions regardless of the 
binning cases. Thus, the oscillation seems to be caused 
by some cosmological origin, not the astrophysical one. 



VI. SUMMARY 

In this work we have re-examined the temperature 
angular power spectrum of the cosmic microwave back- 
ground (CMB) at three frequency bands and for four 
restricted sky directions, in order to explore the origin of 
the fine structures in the power at the multipole range of 
100 < £ < 150. 

We prepared 3,000 mock sky maps of CMB tempera- 
ture fluctuations with anisotropic instrumental WMAP 
noise, and verified the noise and cosmic variance effect 
on the fine structures of the WMAP data. By comparing 
the angular power spectra from each mock data with the 
real data by the WMAP, we found that the probability 
of the fine structures to be realized by chance is about 2 - 
3cr. We checked whether the fine structures of the power 
spectrum depend on the Q, U, and V band maps, and 
found no evidence for the frequency dependence. 

In contrast, we obtained some interesting suggestions 
from the angular power spectra derived from the partial 
skies. We found that the characteristic bump around 
£ = 130 - 140 seen in the all sky angular power spectrum 
is solely attributed to the anomalous power at the South 
East area. We have already known the existence of the 
cold spot as one of the anisotropy structures in South 
East area [111, [HI] ■ We have doubt if the cold spot will 
affect at £ = 140 as a substrucure, thought the cold spot 
is few degree scale. However, we confirm that this effect 
is too tiny to generate characteristic structures at £ = 
140. This result may indicate the existence of unknown 
peculiar structure in that area as the origin of the fine 
structure. The fine oscillating structures found around 
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FIG. 7. Left: Differences in the binned spectra between the predicted angular power spectrum from the ACDM model and the 
WMAP power spectrum with apodization function (9 max = 140° . The black line represents the average value of the 3,000 mock 
data. The deviation from zero manifests the x 2 -like distribution of the angular power spectrum. The boxes are the variances 
estimated from the mock data, as indicated in the figure. Right: The probability contours assuming Gaussian distribution for 
each binned power with no correlation and the WMAP data, as in Fig. 2. 
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FIG. 8. Upper panels: The same as the left panels in Fig. [7] but for the partial skies with binning case I. Lower panels: The 
histograms of differences between the power spectrum of ACDM and those from the 3,000 mock data at I = 138. The red line 
corresponds to the data point of WMAP. The peculiarity of South East sky is clearly shown. 
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FIG. 9. Same as Fig. [5] but for the binning case It and £ — 139 (lower panel). 
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FIG. 10. Same as Fig. [5] but for the binning case HI and I = 137 (lower panel). 
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i = 100 - 120, on the other hand, have no significant 
directional dependences, suggesting that the oscillations 
come from some cosmological origin. 



If the observed feature is not a statistical fluctuation 
but has a primordial origin, a straight but important test 
is to look into the polarization of the CMB and/or the 
distribution of the large scale structure, because they 
should show the same characteristic structure in their 
power spectra. The coming data by PLANCK [33| and 
future galaxy survey such as the Large Synoptic Survey 
Telescope |34[ can be used for the cross check and to 
accept (or reject) the feature with high significance. 
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